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Abstract 

We use the correspondence between the f{R) theory and an Einstein- 
scalar field system to study late-time dynamics of solutions of f{R) theory. 
We discuss how reasonable assumptions on the potential of the scalar 
field lead to restrictions on the function f{R) and use known results for 
the scalar field system to gain results on solutions of the f{R) theory. 
In particular, we prove accelerated expansion at late times for several 
different categories of functions f{R) that satisfy certain restrictions. 

1 Introduction 

Several years ago it was discovered, from the studies of the supernovae data, 
that the Universe is undergoing accelerated expansion and since this time, there 
has been a renewed interest in modified gravity models. There is now such a 
vast amount of evidence for the acceleration of the Universe (supernovae data 
[1], large scale structure formation [2], cosmic microwave background radiation 
[3], weak lensing [4], etc.) that this has led to a great deal of research into 
explanations for the driving force behind this accelerated expansion. Some 
of the simpler theories involve introducing a cosmological constant or some 
kind of dark energy, perhaps driven by a scalar field. However, it may be 
possible to explain the accelerated expansion by modifying Einstein's theory of 
relativity itself. These theories are known as "modified gravity" or "higher-order 
gravity" theories and the particular case we are concerned with here is that of 
"f{R) gravity", where we replace the scalar curvature, R, in the Lagrangian of 
Einstein's theory of relativity with an arbitrary function, f{R), of R. 

The idea of using a non-linear R is by no means a modern one. It became of 
interest in the study of quantum gravity in the 60's and has greatly increased in 
popularity in the last 10 years. There have been many recent papers discussing 
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the viability of /(i?) theories with respect to solar system tests and cosniological 
constraints (for a sample, see [5-13] and refcrences therein). 

Although the literature on f{R) theories is vast, the portion of papers that 
study /(i?) theory from a mathematical viewpoint is much smaller and, although 
it is believed that /(i?) theories may explain the observed accelerated expansion 
of the universe, there are few rigorous mathematical proofs that this is indeed 
the case. This paper will provide some such rigorous proofs for certain cases of 
,f{R) theory. We will not consider a particular form of the function ,f{R), but 
rather prove results on the late-time dynamics of f{R) theories under certain 
restrictions on the function f{R). 

Since the field equations of the f{R) theory are 4th order partial differential 
equations, this makes the theory analytically complicated. However, the theory 
can be related, via a conformal transformation, to Einstein's theory of relativity 
where the matter is governed by a scalar field. Following a suggestion given in 
[14], we can exploit this relationship to gain results on the dynamics of solutions 
to the /(i?) theory using known results on the dynamics of solutions to the 
Einstein-scalar field system. 

The layout of the paper is as follows: In section 2, we give an introduction 
to and derive the field equations for the f{R) theory. Section 3 will deal with 
the vacuum f{R) case. We will first recap the relationship between the f{R) 
theory and the Einstein-scalar field system. We will then consider two different 
restrictions on the potential of the scalar field, for which the late-time dynamics 
of the solutions are known [15,16]. We will translate these restrictions into 
restrictions on the function f{R) and prove results on the late-time dynamics of 
the f{R) theory. In section 4, we will consider f{R) theory with ordinary matter 
satisfying the strong and dominant energy conditions. We will first describe a 
relationship between this /(i?)-matter system and the coupled Einstein-scalar 
field-matter system studied in [17]. We will then use results in [17] on the late- 
time dynamics of the Einstein-scalar field-matter system where the potential 
is restricted in the same way as in the vacuum case to prove results on the 
late-time dynamics of the /(i?)-matter system. 

We will be dealing with cosmological models that are homogeneous and 
isotropic. The underlying spacetimes are the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models. If we assume that the metric of the slices of constant 
time is flat, then the metric takes the form 

ds'^ = -dt'^ + a'^{t){dx^ + dy^ + dz^) 

where a{t) is the scale factor. The scale factor is important for proving accel- 
erated expansion, since it gives a measurement of distances between galaxies. 
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Therefore d > means that the universe is expanding, while a > gives ac- 
celerated expansion of the universe. Another useful parameter is the Hubble 
parameter, H , which is defined by 



We assume c = G = 1. 



2 f{R) theory 



It is well known that the Einstein field equations which describe Einstein's 
Theory of Relativity can be derived from variational principles (see, for example 
[18]). To do this, one considers the Einstein- Hilbcrt action 



C = j Ry/^d^X (1) 

where R is the scalar curvature and g is the determinant of the metric. Taking 
the variation of this action with respect to the metric, ^q^, yields the vacuum 
Einstein equations 

Ra0 - ^Rgo^p = (2) 

where Raf3 is the Ricci tensor. 

We can also add matter to the Einstein equations by choosing a Lagrangian 
of the form 

C = j [R + SttL'^'] ^d^x (3) 

where L^^ is the Lagrangian density for the matter field. Then the Einstein 
equations are given by 

- ^Rgap = SttT^' (4) 



2 - -.y<-^f^ — - a/3 
1/3 



where the energy-momentum tensor, T*^ is defined by 

The matter field must also satisfy equations of motion, which are found by 
taking the variational derivative of the action ^ with respect to the matter 
field. In their general form they arc given by 

where (j) is the matter field. There are many possible choices for the matter 
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and therefore for the energy-momentum tensor. Some describe real, physical 
matter, while others are mathematical tools. For more information on matter 
models see, for example, [19]. 

As previously stated, experimental results over the last ten years have shown 
that the universe is undergoing accelerated expansion. There are several ways of 
adapting Einstein's theory of relativity to account for this accelerated expansion. 
The simplest way is to introduce a cosmological constant to the left-hand-side 
of the Einstein equations. Another, more sophisticated way is to include matter 
in the form of a scalar field. 

However, in this paper, we are interested in an alternative way of explaining 
this accelerated expansion. We will consider the idea of modifying gravity to 
build more freedom into the theory itself. In the study of f{R) theories, we 
replace the R in the Einstein-Hilbert action with a general function of R, to get 
the action 



where a dash denotes differentiation with respect to R, Vq is covariant differ- 
entiation with respect to the metric gap and □ = V^V" is the d'Alembertian 
operator. 

3 The f{R) theory in vacuum 

3.1 fiR) and the Einstein-scalar field system 

The f{R) field equations are 4th order in the metric, compared to the 2nd order 
Einstein equations, which makes them harder to study. Fortunately, we can 
use a conformal transformation to transform them into a system of equations 
which are 2nd order in the metric and then gain results from this simpler set 
of equations, which turn out to be Einstein equations with a specific type of 
matter. 

A reference for this work is given in [20], although it should be noted that 

the signature used in the reference is (H ), compared to the ( h ++) 

signature we use here. 

The conformal transformation that we use is 




(7) 



This then gives rise to the vacuum f{R) field equations 



f'{R)Rap - \f{R)gc.p - VaV^ifiR)) + D{f\R))ga0 = (8) 



gi-^ g where ga/3 = Fga^ 



(9) 
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and where we choose the conformal factor, F, to be the derivative of the function 
/(-R). We find, using the transformation and equation ([8]), that for the g system 
the fohowing equation holds 

-^F~^V^FV^Fg^p (10) 

where Raf3 and R are the Ricci tensor and scalar curvature built from the metric 
g and where Vq is covariant differentiation with respect to g. If we then make 
the choices 

and 

VW-^ — ^^^^ (12) 
we find that this equation becomes 

Rcp - ^Rgap = Stt (^V„0V;30 - iv^0V''<?ig„/5 - y(0)5a/3^ (13) 

The right hand side is precisely the form for the energy-momentum tensor of a 
scalar field, (f>, with potential V((j)) and so we find that the conformal metric, g, 
is a solution to the Einstein equations with matter described by a scalar field. 
If we take the trace of the f{R) field equations (O, we get 

f{R)R-2f{R) + 3a{f'{R)) = 

or, since F — f'{R), 

FR - 2f{R) + iUF = 
Multiplying through by j^y^^F^"^ gives 



1 Wtt 
16^V ~3~ 



FR - 2f{R) 
F^ 



Recalling the definition of the potential, V, of the scalar field shows that the 
derivative, V {(/)), with respect to (j) is exactly the 1st term on the left-hand-side. 
We can then make the conformal transformation g i— > g as before. We get 

OF = FDF -V'^FVc.F 
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Since F and 4> are related hy F = ev 3 "f", this then gives 



^^F^Uc^ 
3 ^ 

We can then substitute this information into equation p4p above to show that 

V\<j,) = 00 (15) 

which is exactly the equation of motion for the scalar field, <f>. 

From the definition of the scalar field in terms of the function f'{R) 



(16) 

it is not clear whether we can express the scalar curvature, i?, in terms of the 
scalar field, (j), for all choices of f{R). We would like to be able to do this in 
order to express the function of i? 

1 fiR)-Rf'iR) 

in terms of <f>, since V represents the potential of the scalar field. We also require 
that V is smooth. In order to achieve this, consider as a map : R R given 

by 

In order to ensure that R can be defined in terms of (/>, we require that the inverse 
of this map exists and is smooth. By the inverse function theorem, we see that 
this is the case if (/)'(i?) 7^ for all R. So we must insist that f"{R) 7^ for all 
R. In other words f'{R) is one-to-one and therefore either strictly increasing or 
strictly decreasing. 

Hence there is a nice correspondence between the vacuum /(i?) theory and 
the Einstein-scalar field system, where solutions of these two different theories 
are related by a conformal factor, F , which is the derivative of the function 
f{R) which appears in the f{R) theory. 

Therefore we can use known results on the dynamics of solutions of the 
Einstein-scalar field system to find results on the dynamics of solutions of the 
f{R) theory. For the Einstein-scalar field system, the dynamics of solutions 
depend on the form that the potential of the scalar field takes. There are several 
known results for this system in the case of a homogeneous space-time, where 
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different assumptions on the potential have been made. In order to relate this 
to the f{R) theory, we must translate these assumptions on the potential into 
assumptions on the function f{R). We will then gain results on the dynamics 
of the f{R) theory which depend on the function f{R). The first case we will 
consider is the case where the potential admits a strictly positive minimum. 

Since we will be interested in accelerated expansion, it is important to know 
the relationship between the Hubble parameter of the f{R) system and that of 
the Einstein-scalar field system. From the relationship between the metrics, we 
find that the scale factors are related by 

a = Fa (18) 

and so the Hubble parameters are related by 

H=~^ = ^:^^U^^H+^ (19) 
a Fa a F F ^ ^ 

3.2 The potential admits a strictly positive minimum 

Recall the relationship between the scalar field and the function f{R) given 
in We see from this that we must have f'{R) > 0. The potential F(0) of 

the scalar field and the function f{R) are related by equation ([T7|). 

We would like to consider cases of the scalar field where the potential is 
strictly positive and has a strictly positive minimum, since these cases lead to 
accelerated expansion of the universe which does not stop [21], and to study 
how the insistence of a strictly positive minimum gives rise to restrictions on 
R, f{R) and the derivatives of f{R)- 

To describe the asymptotic (latc-time) behaviour of a solution of the f{R) 
theory which satisfies these conditions, we follow the paper in [15], which gives 
results on the asymptotic behaviour of a solution to the Einstein-scalar field 
system, and apply these results to the equivalent f{R) theory. 

Consider first the case where the potential, V, is strictly positive. This 
condition becomes, in terms of R and /(i?), 

f{R) - Rf'iR) < (20) 

for all R. 

We can already prove some results on the f{R) theory with this restriction 
in place. 

We first prove an analogous theorem to Theorem 1 in [15] for the scalar field 
system, which will be more useful for the translation to f{R) theory. 



7 



Theorem 1. Consider a solution, g, of the Einstein equations of Bianchi type I- 
VIII coupled to a nonlinear scalar field with potential V of class defined on the 
(possibly infinite) interval (a, (3) and satisfying the following two assumptions: 

1- yi4>) ^ ^0 for a constant Vq > 

2. V tends to oo as 4> tends to a or (3 
If the solution is initially expanding (H > Q) and exists globally to the future, 
then fort — + oo, the quantities R andi/^— (87r/3)[0^+V^(^)] decay exponentially. 
The potential V{(j)) converges to some positive constant Vi, V'{(f>) and 
H (SttVi/S)^^^, where H is the Hubble parameter for the Einstein-scalar 
field system. 

Note that this theorem will also hold in the presence of other matter satis- 
fying the dominant and strong energy conditions, as in Theorem 1 in [15], but 
this is not needed for the purpose of this paper. 

Proof. Assumption 2 of the theorem above {V ^ oo as (j) ^ a or P), together 
with the information that V is bounded to the future (by conservation of energy) , 
tells us that cj) must remain bounded away from a and P as t oo. Therefore, 
since V is smooth, the derivative of V with respect to cj) must also be bounded 
to the future. This information is sufficient to replace the second assumption 
in Theorem 1 of [15], since the proof only requires that V is bounded to the 
future. The third assumption in Theorem 1 of [15] is required to show that, if 
(f) — > ±oo as t oo, we know that V' converges. Now, since V is defined on 
(a,/3), then may only converge to ±oo as < — *■ oo if one of a or /3 is infinite. 
However, we know from above that (j) is bounded away from a and /3 as t ^ oo 
and so (j) must not converge to ±oo as t oo. Therefore, this assumption is 
no longer required. With these alterations, the proof of Theorem 1 above now 
follows directly from the proof of Theorem 1 in [15]. □ 

To translate this theorem to the f{R) theory, we therefore need some condi- 
tions on the function f{R) which will ensure that the potential blows up at its 
endpoints. 

We know already that the function f'{R) is strictly positive and is one-to- 
one. For the purpose of this paper, we will restrict ourselves to the case where 
f'{R) is defined on the whole real line. The situation where f'{R) has a smaller 
domain of definition can be treated in a similar way. We assume for now that 
f'{R) is strictly increasing. Then, considering the graphs of f'{R) against R 
and (j) against f'{R), we find that both are strictly increasing. Suppose that 
the range of f'{R) is (a, 5) (the biggest range posssible being (0,oo)). Then (j) 
can take values between \ / tI- In a and a / -S- In b. So the domain of definition of 

the potential V is iJ In a, J In b J (the biggest possible domain being the 
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whole real line) . We are interested in how V behaves as cj) approaches y In a 

and In b and we sec that this is equivalent to studying the behaviour of V 

as R tends to — cx) and oo. 

So we must consider how V behaves as i? ±00 and show what assumptions 
must be made in order that V blows up. This will depend on the range of ,f'{R)- 

Consider the case where f'{R) tends to a (positive) constant a as i? ^ —00. 
Then we can write, using the Landau notation [22], f'{R) = a+o(l) as i? ^ —00. 
Integrating gives f{R) = aR + o{R) as i? — > —00 and multiplying f'{R) by R 
gives Rf'iR) ^ aR + o{R). Subtracting Rf'{R) from f{R) shows that the 
leading order terms cancel and we get f{R) — Rf'{R) = o{R). Recalling the 
definition of the potential in terms of the scalar curvature R, we find that as 
R —00, 

167r(a + o(l))2 

This may or may not blow up, depending on the exact form of f{R) — Rf'{R) as 
R —>■ —00. In order to ensure that the potential tends to infinity as i? ^ —00, 
we require that f{R) — Rf'{R) —00 as R^ —00. 

Now consider the case where f'{R) ^ as i? ^ —00. In Landau notation, 
this means that f'{R) = o(l) as i? ^ —00 and so f{R) — o{R) as i? ^ —00. 
So we find that Rf'{R) ~ o{R) as i? ^ —00 also. Since we know that f'{R) 
tends to zero as R tends to negative infinity, we see that the denominator of 
the potential tends to zero. Therefore, provided the numerator does not tend 
to zero faster than the denominator, we will find that the potential blows up 
as R tends to negative infinity. One way to ensure that this is the case is to 
restrict the function f{R) so that the expression f{R) — Rf'{R) is bounded away 
from zero. (It should be noted that this is not a necessary condition, but is a 
sufficient one.) 

This covers the situation as R — *■ —00. Now let us consider what happens 
as i? — * 00. 

First, suppose that the endpoint of the range of f'{R) is finite and equal to 
b. (So f'{R) ^ 5 as i? ^ 00). Then f'{R) =b + o(l) for R ^ 00. Integrating 
gives f{R) = bR + o{R) and we also have Rf'{R) = bR + o{R) and {f'{R)f = 
{b + o{\)Y as R ^ 00. Therefore we have the same situation as in the case for 
f'{R) — > a as i? ^ —00. Once again, we find that the potential will blow up 
provided f{R) — Rf'{K) —00 as R ^ 00. 

Now suppose that f'{R) tends to infinity as R tends to infinity. The situation 
is now more complicated. If ,f'{R) —>■ 00, then f{R) must also tend to infinity. 
So we have the situation that R, f{R) and ,f'{R) all tend to infinity. We know 
that the denominator of the potential tends to infinity as i? — > 00. So we must 
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insist that the numerator tends to negative infinity as R oo. However, this 
is not enough to ensure that the potential blows up. We must also insist that 
the numerator tends to negative infinity faster than the denominator tends to 
infinity. This requires that (/'(i?))^ = o{f{R) - Rf'{R)). 

The case where f'{R) is strictly decreasing can be treated in a similar way 
and we find the following restrictions on R, f{R) and f'{R)- 

(Fl) If f'{R) ^ a as i? ^ oo or — oo, we require f{R) — Rf'{R) — oo as 
i? ^ oo or — oo. 

(F2) If /'(i?) as i? oo or -oo, we require that f{R) - Rf'{R) is 
bounded away from zero as i? ^ oo or — oo 

(F3) If/'(i?) ^ oo as i? -> oo or -oo, we require that/ -oo 
as i? oo and that {f'{R)f = o{Rf'{R) - f{R)). 

We are now ready to state the theorem for the f{R) theory. 

Theorem 2. Consider a solution, g, to the f{R) field equations in a Bianchi 
type I-VIII space-time, where f{R) satisfies the following properties: 

1. f{R) - Rf'{R) < for all f'{R) > 

2. f{R) - Rf'{R) satisfies the three restrictions (Fl), (F2), (F3) 

If the solution satisfies H + ^ > initially and exists globally to the future, 
then for t —^ oo, the quantity — tends to some positive constant 

K, '^^^^fi(ji))2'^^ tends to zero and H tends to {^)^^^, where H is the Hubble 
parameter. 

Proof. From the relationships between the function f{R) of the f{R) theory 
and the scalar field (j> of the Einstein-scalar field system, between f{R) and the 
potential, V{(p), and between the Hubble parameters of the two systems, we see 
from the discussion above that, if g satisfies the assumptions of this theorem, 
then g ~ Fg is a solution of the Einstein-scalar field equations satisfying the 
assumptions of Theorem 1 above. Therefore, we can use the results of the latter 
to show that, if the solution is initially expanding {H > 0) and exists globally 
to the future, then for t — > oo, the potential V{(j)) converges to some constant 
Vi, V'{(j)) — > and the Hubble parameter, H, of the Einstein-scalar field system 
tends to (^)'^'. 

Translating this back to the f{R) system, using the relationships between 
f{R) and (p and between f{R) and V{(f>), gives the first two results of Theorem 
2 above, where the constant K is equal to IGttVi (and is therefore positive, since 
Vi is positive). Recall that the Hubble parameters are related by 

F 

H = H + - 
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We can use the information from the proof of Theorem 1, for which we refer 
the reader to the proof of Theorem 1 in [15], that ^ as t oo. From the 
definition of in terms of F (equation (jlip). we get 



and so we find that |i ^ as t — > oo and hence that the Hubble parameter H 
is equal to H at late times. This concludes the proof of Theorem 2. □ 

The theorem therefore states that, as < ^ oo, the quantity '^^'^^f~(B)Y'^'^ 
tends to zero and the quantity ""^^^p^^^r^ tends to a constant K. Suppose 
now that any points where the quantity '^^'' ffi{B)Y^^^ equals zero are isolated. 
Then this, together with the information that this quantity tends to zero as 
< — > oo, tells us that R must converge as i ^ oo. But — ^ ^^(f~(ji))'i^^ — > oo as 
R ±oo (from the assumptions (Fl), (F2), (F3) of Theorem 2), together with 
— '^^^/'(fl()2^^ — !• as t OO, shows that R must converge to a finite limit, Rc 
say, as t — > oo. So we find that a.s t oo, R ^ Rc with 2/(i?c) — Rcf'{Rc) ~ 0. 

Considering again the equation for V in terms of f{R), given in (fT7|) and 
taking its derivative with respect to 0, tells us that 



dV _ 1 /IGtt 
'M ^ 167rV ~3~ 



2f{R) - Rf'jR) 



(22) 



We therefore see that the limit, i?c, of i? as t oo corresponds to a critical point 
of the potential, V . Furthermore we can use the fact that 2/(i?c) = Rf'{Rc) to 
find the value of V at this point, and therefore the constant K of Theorem 2, 
in terms of Rc ■ We find that 

K= 



4/(i?c) 



which shows that H tends to a positive limit and so we have accelerated expan- 
sion at late times for the solution of the f{R) theory. 

To gain more detailed information about the late-time asymptotics of the 
solution, we impose further conditions on the potential, V. We assume now that 
V tends to a strictly positive non-degenerate minimum as t oo. To discover 
how this restricts the function f{R), we use (fTS)) and (fT7|) to find the second 
derivative of the potential with respect to 4>. We get 



d^V _ 1 
Id^ ~ 3 



if'iR)r + RnR)f"{R) - ^f{R)f"{R) 

{fiRWinR)) 



(23) 
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where f"{R) denotes the second derivative of / with respect to R. The point 
Rc corresponds to the critical point of the potcntiaL Recalhng that this point 
must satisfy 2/(i?c) — Rcf'{Rc) and substituting this into wc get 



d^V _ 1 
del? ^ 3 



1 Rc 



f"{Rc) f'{Rc) 



Then, in order that the critical point is a non- degenerate minimum, we require 
that the second derivative is positive. This requires that the term in brackets is 
positive. Since the expression f{Rc) — Rcf'{Rc) must be negative (from (PO]) ) 
and we have the equality 2/(i?c) = Rf'{Rc), we find that Rc and f{Rc) must 
both be positive. Then in order that the term in brackets is positive, we need 

f'iRc) > 

and 

f'iRc) < ^'^""'^ 



Rc 

So, in order that the critical point is a non-degenerate minimum, we must have 
that the second derivative of / with respect to R at this point lies in the interval 

■f'"^'' 



We can now use Theorem 2 in [15] to get more detailed information on the 
asymptotics in the future. 

Theorem 3. Let g be a solution to the f{R) field equations satisfying the hy- 
potheses of Theorem 2. 

Then, if Rc is the limit of R as t oo, and 

< f'iRc) < 

Uc 

we find that the quantities f'iRc), ~ ~ ^ '^'"^'^ ^ ^ decay ex- 

ponentially as t oo, where Hi = \J^- We also find that /'(-Rc)e~^^^*gQ/3 
converges to a limit. 

Proof. The assumption made on the second derivative of /(i?) as i ^ oo is 
equivalent to the assumption made on the second derivative of the potential in 
Theorem 2 in [15]. It is then seen that the conformally transformed metric, 
5q/3 = Fgafs, will satisfy all the assumptions required for that theorem and 
hence that the conclusions will follow for the metric (ja/s . Using the information 
gained in Theorem 2 of this paper and the discussion after, we see that the 
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results of Theorem 2 in [15] can be transferred onto the metric, g, of the f{R) 
theory to give the conclusions of Theorem 3 above. □ 



Remark 1. It should perhaps be noted here that the condition that the minimum 
of the potential is strictly positive is very important for the accelerated expansion 
result. To see this, consider a function of the form 

f{R) ^R + aR^ 

which gives a potential 

487ra V / 

that admits a zero minimum at (f) = 0. The paper [23] has studied fiat and open 
FLEW models with a scalar field that has a potential with a local zero minimum. 
It is shown that the solution is ever expanding, but that and H asymptotically 
approach zero and therefore the solution does not undergo accelerated expansion 
at late times. Therefore, as discussed in section 4 of [24], the corresponding 
f{R) function which gives rise to an Einstein-scalar field solution with this type 
of potential cannot provide a mechanism for late accelerating expansion of the 
universe. 

To finish, I would like to bring attention to an example [10] in the literature 
of a function f{R) that fits the criteria of this section and is therefore a solution 
of the f{R) theory that undergoes accelerated expansion. 

In the paper, the function studied is 

/(i?) ^ ii? + J_ log (cosh(ai?) - tanh(6) sinh(ai?)) (24) 

where restrictions are placed on the parameters a and b by requiring a > and 
b > 1.2. For the purpose of this paper, we assume a > 0, but will not place 
restrictions on 6 until they are needed. The derivative of f{R) is given by 

/'(i?) = i [l + tanh(ai?-6)] (25) 

and this clearly satisfies the condition that f'{R) > for all R. Note that the 
range of f'{R) is (0, 1) and hence, from (|16|). that is only defined for (— oo, 0). 
This means that the domain of definition for the potential is (— oo, 0). 

The condition that the potential is strictly positive, given by equation (|20p . 
becomes 

^ log (cosh(ai?) - tanh(6) sinh(ai?)) - i?tanh(ai? - 5) < 
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The expression on the left-hand-side tends to —oo as i? — > ±00 and contains 
only one turning point, a maximum, which takes the value zero at i? = 0. So 
we find that the inequality is true for all R, except the point R = 0. 

Now, the function f'{R) behaves in the following way as i? — > ±00. As 
R —00, f'{R) and as i? ^ 00, f'{R) — > 1. So in order that the second 
assumption in Theorem 2 (corresponding to the blow-up of the potential at its 
endpoints) holds, we require that f{R) — Rf'{R) —00 as R ±00. But we 
have already shown that this is the case, and so the function satisfies the second 
condition of Theorem 2. 

The fact that we do not have that f{R) — Rf'{R) < for all R seems to cause 
some problems. However, we can get round this by considering the conservation 
of energy, which tells us that if a solution falls into a strictly positive local 
minimum of the potential, it will get trapped there. If this is the case then 
we are free to change the shape of the potential away from this local minimum 
without affecting the solution. Hence we need not worry about the behaviour 
of the potential away from the strictly positive local minimum. 

The discussion following Theorem 2 then tells us that, as, t ^ 00, R Rc 
with 2/(i?c) = Rf'{Rc)- We can apply this to the function under consideration 
to find the possibilities for the limit, i?^; of i? as t ^ 00. So let us study the 
zeroes of the function Q{R) = 2/(i?) - Rf'{R). We get 

QiR) = f + - log (^^^^^^^) - f tanh(ai? - b) (26) 
2 a \ cosh(o) / 2 

Consider what happens to Q{R) as R tends to plus and minus infinity. As 
R -00, cosh(ai? -b)- ^e"'-''"-'''^ and tanh(ai? - b) ^ -1, so 

Q(i?) ^ - i log 2 + - - i log(cosh(&)) 
a a a 

This expression is negative for any choice of 6, since a > 0. To see this, consider 
the graph of this expression against b. We see that it tends to —00 as 6 tends 
to —00 and zero as b tends to 00 and has no turning points. Therefore, as 
R tends to —00, Q{R) tends to a negative limit. As R tends to plus infinity, 
cosh(ai? — 6) — ^e°-^^^ and tanh(ai? — &) ^ 1, so we find that Q{R) tends 
to infinity. 

Now, at i? = we see that Q = 0, which gives one possible choice for Re- We 
can also see that the function Q is increasing between minus infinity and zero, 
so Q does not have any zeroes for i? < 0. We would like to know if there are 
any zeroes for positive R. To investigate further, consider the second derivative 
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of Q, which is given by 



Q"{R) tanh(ai? — b) 

cosh^ (ai? — b) 

For & < 0, the second derivative is always positive and so the function Q{R) is 
convex and wiU clearly have no more zeroes. In this case the only possibility 
for Rc is Rc — 0. But if this is the case then we find that the constant K in 
Theorem 2 is also zero and this docs not give rise to a solution which undergoes 
accelerated expansion. So if we want a solution that gives accelerated expansion, 
we require more options for Rc and hence a zero of Q{R) at some i? > 0, where 
the constant K of Theorem 2 will then be positive. Therefore we must assume 
& > 0. Clearly, at i? = ^, we have 

Q"{R) = 

This tells us that the graph changes shape at ^. For R < ^ the second derivative 
is negative and for i? > | it is positive. Now, at i? = ^ , the function Q is given 

by 

Q(.R) = ^{b^2logicoshm 
la 

Clearly if the value of Q(R) is negative at i? = ^, then there must be 2 more 
zeroes of the function Q (since at i? = 0, the function Q is zero and increasing 
and as i? oo, Q ~* oo). Since a > 0, we must investigate for what values of b 
we have 

6- 21og(cosh(6)) < 

It can be shown that choosing b > 2 log 2 ensures that this inequality holds. 
To see this, note that if 

5- 21og(cosh(6)) < 

then 

gb ^ glog(cosh2(6)) 

and so 

2ef - - 1< 
Making the substitution x ~ e^/^ gives 

2e¥ - e^** - 1 = 2a;^ - - 1 

the right-hand-side of which can be factorised as 

2x3 - a;^ - 1 = (x - i)[-x^ + + x + 1) 
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Since we have assumed that 6 > we see that x > 1 and so x — 1 > hence we 
see that 

6-2 log(cosh(&)) < 4=^ + x^ + x + l<0 

and smce x ~ e**/^ we must have positive x. Studying the polynomial —x^ + 
x'^ + X + 1 shows that its positive root lies between a; = 1 where the polynomial 
is positive and x — 2 where the polynomial is negative. Hence, if we take 
X > 2, then we are sure that the polynomial is negative. This corresponds to 
b > 2 log 2. 

Therefore, if 5 > 2 log 2, then the function Q{R) has 2 further zeroes and so 
we have another 2 possibilities for the limit, Rc, of i? as i ^ oo (since we have 
already ruled out the case Rc = 0). 

For Theorem 3, we require an extra condition on the f{R) function, that 

o<f"{Rc)<^^ 

We see that this is equivalent to Q'(R) > 0, since from the definition of Q{R) 
we get 

Q'{R)^r{R)-Rr(R) 

Using all the information we have gained about Q{R), we can see that the 
third zero of Q must have positive gradient and the second must have negative 
gradient. (To see this, consider the sign of Q'{R) immediately before and after 
i? = ^). Therefore, if the limit of i? as t oo is the third possibility given, 
we see that this Rc satisfies the conditions of Theorem 3 and so this choice of 
f{R) admits a solution of the f{R) theory which gives rise to the accelerated 
expansion described in Theorem 3. 

To give a concrete example, we choose, as in the plots given in [10], a = 2 
and 6 = 1.5. If we consider the graph of the function 2f{R) — Rf'{R), we 
see that 2f{R) - Rf'{R) has 3 zeroes at i? = 0, i? ~ 0.4818758472 and at 
R ~ 1.385168809. So we have 3 candidates for the choice of Rc- The first 
of these gives K = in Theorem 2 and so does not give rise to accelerated 
expansion. In the second and third cases, we get 

f"{R) = ^ 



cosh^(2i?- 1.5) 
and 

^ = ^(l + tanh(2i?-1.5)) 
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Wc find that the second case gives 



(which does not satisfy the assumptions of Theorem 3) but the third case gives 

f'{R) 



.f"{R) < 



R 



Therefore this case satisfies the assumptions of Theorem 3, as required. See the 
paper [10] for a plot of the potential. 



3.3 Intermediate inflation 

We have considered the case where the potential is strictly positive and has a 
minimum. The next step is to allow the potential to tend to zero as the scalar 
field tends to infinity. In particular, we would like to consider potentials which 
satisfy the following two conditions: 

1. > with V{(t)) as -> cx) 

2. V'{(l)) < 

The paper [16] gives some results on the dynamics of scalar fields which satisfy 
these conditions. It is known that accelerated expansion will occur provided the 
potential does not decay too rapidly. Therefore it is necessary to restrict the 
potential further, e.g. placing a restriction on the ratio V /V will ensure accel- 
erated expansion at late times. The paper [16] places the following restriction 
on V'/V: 

3. V'{(t))/V{(j)) ^ as cx) 

Then Theorem 1 of [16] shows that the solution undergoes accelerated expansion 
at late times. We can use this theorem to prove results on the late-time dynamics 
of the f{R) theory. 

To proceed, note that the assumptions 1, 2 and 3 above on the potential 
of the scalar field can be translated into assumptions on the function f{R) by 
means of equations (fT7|) and ((22|) . We get 

(fl) /(i?) - Rf{R) < with ^ as fiR) ^ ex. 

(12) 2f{R) - Rf'{R) < for all R 

(f3) ^g) ^ 2 as f'{R) ^ cx). 

We assume (analogously to the paper [16]) that we are in a situation where 
f'{R) ^ oo as t oo. This assumption will be required for the rest of this 
section. We can now use Theorem 1 of [16] to prove the following result on 
solutions of the f{R) theory that satisfy conditions (fl), (12) and (f3). 
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Theorem 4. Consider a spatially flat, homogeneous and isotropic solution of 
the f{R) field equations satisfying conditions (fl), (f2) and (fS) above. Suppose 
that the solution satisfies H > and y > initially, where F = f'{R) and H 
is the Hubble parameter. If the solution exists globally to the future, then the 
scale factor, a, satisfies a > for t sufficiently large and the solution undergoes 
accelerated expansion at late times. 

Proof. Firstly, we have shown that, if g is a sohition to the f{R) field equations, 
given in satisfying conditions (fl),(f2) and (f3), then g ~ Fg is a solution 
to the Einstein-scalar field equations, given in (fT5|) . where the potential of the 
scalar field satisfies conditions 1-3 above. Now, recall that the relationship 
between the Hubble parameter, H, of the Einstein-scalar field system and the 
Hubble parameter, H, of the f{R) theory is given by (fT9|) 

F 

H = H + - 

and that the scalar field and the function F are related by H 



— Ini^ 
IGtt 



Therefore cj) = y j; and we find that the assumptions in the theorem above 

that H > and -p > initially, tell us that H > initially and that (j) > 
initially, which ensures that all assumptions of Theorem 1 in [16] are satisfied. 
Therefore the results of that theorem hold for the solution, g, of the Einstein- 
scalar field system. In particular, a > for sufficiently large t and the solution 
undergoes accelerated expansion. 

We would now like to show that the scale factor, a, of the f{R) theory also 
satisfies a > for sufficiently large t. To do this, we note that this is equivalent 
to proving H + > 0. We know from Theorem 1 of [16] that the inequality 
H + > (which is equivalent to |^ < 3/2) holds for sufficiently large 
t. So we consider the relationship between H and H. We know already that 
H = H + j: and that (j) = \J j^T^- During the proof of Theorem 1 in [16], it is 
noted that 

3i/2 _ 3i/2 



So we find that 



— > as > 1 and > 1 as i ^ cxd 



as t — > cx) 



H 

But from the first assumption on the potential, we know that \^ ^ as t — > oo 
and so, since |^ ^ 1 as t ^ oo, this means that H ^ a.s t ^ oo. Hence, 
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from above, we must also have that — > as t ^ oo. 

Therefore wc find that ^ ^ as t ^ oo. So, from the relationship between 
H and H, we see that H behaves like H at late times. Now consider H. This 
is related to the H of the Einstein-scalar field system by 



H 



We know already that the third term tends to zero as t tends to infinity. From 
the relationship between F and (f> we find that 




F _ IGtt -2 /167r •; 



F 3 " V 3 " 

Once again, we know that (/> ^ as i — *■ oo. We also know, from [16], that the 
following equation holds for the Einstein-scalar field system: 

4>= -3H(f>-V'{(l}) (27) 

and that — > and (/) — > as i — > oo. From the first and third assumptions on 
the potential of the scalar field, we see that V ^ and V /V ^ as t oo, 
which tells us that V^' as t oo. So we find that (j) must also tend to zero 
as t tends to infinity. Thus we find that 

Therefore the behaviour of ij + as t ^ oo mimics that oi H + and, 
in particular, H + > at late times. This concludes the proof. 

□ 

We would now like to give an example of a function f{R) that satisfies the 
conditions (fl), (12) and (f3) above. To find such an example, we note that the 
function 

f{R) = i?2 

gives rise to a scalar field with constant potential, whereas the function 

f{R) = i?" 

with n > 2 gives a scalar field with exponential potential. We are looking for 
something that lies in between these two functions. Consider 

f{R)^RHogR 
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defined on (e~i/2,oo) so that f'{R) > 0. Then f{R) > with f'{R) ^ oo as 
R ^ oo. Now for the first condition, (fl), we get 

f{R) - Rf'{R) = -R^ logi? - R^ 

which is negative for all R in the domain and 



f{R) - Rf'iR) 



-1 - 



log/?. 



as f'{R) -> oo 



if'iRW (2(logi?)l/2 + (logi?)-l/2)2 

For the second condition, (f2), 

2f{R) - RfiR) = -R^ 
which is negative for all R in the domain and for the third condition, (f3), 

i, J =2+ >2asR->oo 

fiR) logi? 

as required. 

Therefore, this particular choice of fiR) satisfies the conditions of Theo- 
rem 4 and is another example of an fiR) theory which undergoes accelerated 
expansion at late-times. 

In the paper [16], it was shown that the third condition on the potential that 
V' /V ^ as ^ oo can be relaxed to 

limsup(-y7V^) < 4^^F 

without altering the result that 5 > 0. This then gives accelerated expansion 
for a more general type of potential. In particular, this includes the case of the 
exponential potential. 

We can similarly relax the condition (f3) on the fiR) theory to prove ac- 
celerated expansion results on a more general type of fiR) function. Recall 
that 



V mn f2fiR)- Rf'iR) 



V V 3 \ fiR) -Rf'iR) 



l6Tr ( fiR) 



3 \fiR)- Rf'iR) 
Then our alternative condition for the fiR) theory is 
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We can now prove the following theorem: 



Theorem 5. Consider a spatially flat homogeneous and isotropic solution to 
the f{R) field equations satisfying conditions (fl) and (J2) above with 



and H is the Hubble parameter. If the solution exists globally to the future, 
then the scale factor, a, satisfies a > for t sufficiently large and the solution 
undergoes accelerated expansion at late times. 

Proof. As in the proof of the previous theorem, we know that if 5 is a solution 
to the f{R) field equations, given in Q, satisfying the conditions of Theorem 
5, then g ~ Fg is a solution to the Einstein-scalar field equations, given in 
([T^ . where the potential of the scalar field satisfies conditions 1 and 2 above 
on the potential, together with the condition limsup(— F'/F) < Ay/n. We also 
know, as in the previous theorem, that, ii H > and > initially, then 
H > and > initially. Therefore, the solution g of the Einstein-scalar field 
system satisfies the assumptions of Theorem 2 in [16] and wc know that, for this 
solution, limsup (^^§y^ < 3/2, which implies H + > 0, which gives a > 0. 
Now recall that 



By the definition of the scalar field in terms of the function F = f'{R), we have 




Suppose that the solution satisfies H > and y 



> initially, where F ~ f'{R) 




and 





In the proof of Theorem 2 in [16], we find that 




as i — *■ 00. Then, since 




we find that 



limsup — < 
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as t — > oo. Now limsup (^|^^ < | as t — ^ oo together with y — > as t ^ oo 
(from the first assumption of Theorem 2 in [16]) imphcs that 

H ^ sls t ^ oo 

Then, this information together with hmsup (^"^^ < as i ^ oo imphes 

that 

^ as t ^ c» 

Hence we know that 

F 

— ^ as i ^ oo 
r 

As before, we know that (p = —iHcj) — V'{4)) and that and as 

t — > oo. We also know, from the assumptions of Theorem 2 in [16] that 

limsup {~V' /V) < A^/tt as i ^ oo 

and that ^ as t — ^ oo. Therefore we must also have that y ^ as i oo. 
Hence, from above, 

— > as t — ^ oo 

In summary, we have shown that |^ — s- and ^ as t oo and hence 
that H + behaves like H + for large t and in particular that a > at 
late times. □ 

There is an important example of the f{R) theory that satisfies the condi- 
tions of Theorem 5, namely 

f{R) = ai?" (28) 

defined on the interval (0, oo), where a is a positive constant and n > 2. Then 
/(i?) - Rf'{R) = ai?"(l - n) is negative for all R G (0, oo) since n > 2. Now 

f{R)-Rf{R} ^ 1-n 2-n 
(/'(i?))2 an^ 

which, again since n > 2, tends to zero as f'{R) tends to infinity. So the 
condition (fl) is satisfied. Condition (f2) is trivially satisified by noting that 
2/(i?) - Rf'{R) = ai?"(2 - n), which is negative for all R e (0, oo). The final 
condition that 

i^"^^'^p(/(i?)^-?m)<^"' 
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is shown to be true by noting that 



fjR) ^ 1 

f{R) - Rf'iR) 1 - n 

which is, in fact, negative since n > 2. Hence this particular form of /(i?) is an 
example of a function which satisfies the conditions of Theorem 5. 

This is an important example, since from the definitions of V and in terms 
of R, we can show that this choice of f{R) corresponds to the following form of 
potential: 

where A — (n^) ' ^^'^ since it can be shown that A < \/2 for all n > 2, we 
see that this is precisely the exponential potential described in [16] and in [25]. 



4 f{R) with matter 

Up to now, we have only considered cases where we have vacuum solutions to 
the /(i?) field equations. We would like to prove some results in the case where 
we add matter. Suppose that g is now a solution to the ,f{R) field equations 
with matter described by the energy- momentum tensor T^p. Then the action 
we consider is given by 

£ = y" [/(i?) + 8TTL'^']^d'^x (30) 

where L^^ is the Lagrangian density for the matter field. Taking the variation 
of this action with respect to the metric gives the field equations for the f{R)- 
matter system: 

f'{R)Rap - IgcpfiR) + □(/'(i?))5a/3 - V„VM/'(i?)) = SnT^f; (31) 
where the energy-momentum tensor, T^^ is defined by 

rpM , ^ tM „ /qoA 

"'^""9^ 2 ^'^'^ ^ ' 

The matter must also satisfy an equation of motion, which will depend on 
the type of matter used. In general, the equation of motion will look like 

^ - V" (^^\ (33) 
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where ip is the matter field. 



4.1 f{R) with matter and the coupled Einstein-scalar field- 
matter system 

In the case of f{R) tlicory witliout matter, wc fomid tliat the f{R) theory has 
a one-to-one equivalence with the Einstein-scalar field system. We can similarly 
show that when we add matter to the f{R) theory, there is a correspondence 
between the /(i?)-matter system and the coupled Einstein- scalar field-matter 
system, which has been studied in [17]. The paper [17] studied the coupled 
Einstein-scalar field-matter system and proved results on the late-time dynamics 
for certain choices of the potential of the scalar field. By use of a conformal 
transformation we show here that there is a correspondence between the f{R)- 
matter system and the coupled Einstein-scalar field matter system and use the 
results in [17] to prove results on the late-time dynamics of solutions of the 
/(i?)-matter system. 

Let g be a solution to the coupled Einstein-scalar field-matter equations. 
Then g must satisfy the following system of equations: 



where Raf} and R are the Rieci tensor and scalar curvature built from the metric 
g, Vq is covariant differentiation with respect to g, (f> is the scalar field with 
potential V{(j)), T^p is the energy-momentum tensor for the ordinary matter 
with T^^ its trace, V'{(f>) is the derivative of V with respect to (p, □ is the 
d'Alembertian operator and C{(j)) and c(^) are the coupling constants. 
The energy-momentum tensor can be written as 



where Saf3 is the spatial component of T^^, ja is the current density and p is 
the energy density. 

We can now show the correspondence between the /(i?)-matter system and 
the coupled Einstein-scalar field-matter system. To do this, we let g be a solution 
to the /(i?)-matter field equations given in (|3ip . The energy- momentum tensor 



Ra/B — X^5ct/3 



87r[V„<^V;30 - -Vy^V''^go.p - y(0)ga/3 (34) 




□0-r(0) 




f^pdx'^dx'^ = S -j®di-di®j + pdt® dt 
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T^p is given (similarly as above) by 

T^jidx^dx'^ := S - j (g) dt - dt (g) j + pdt (g) dt (37) 
We then make a conformal transformation as follows (c.f. section 5 in [17]). 

g^g:= Fg t^t:= f F^^'^dt 

J to 

S :=S j ^ J := F-'^/^j P ^ P ■= F''^ p 

Then we see that the components of the energy- momentum tensor, T^^, in the 
g coordinate system are the same as the components of the energy-momentum 
tensor, T^^, in the g coordinate system. 

Using this transformation and equation (|3ip . we find for the g system that 
the following holds (c.f. equation pO]) for the vacuum case) 

Rcp - \R~gc,p = »7:F-'f^f, + if{R)F-^o.f3 - \F~^Rgo.p 
+ \f-''V C.FV pF - ^F-^V^FV^Fg^0 

where Rafi and R are the Ricci tensor and scalar curvature built from the metric 
g. If we then make the choices (as in the vacuum case) 



IGtt 

and 



^ Ini^ (38) 



as well as the choice 

C{<f) = (40) 
we find that this equation becomes 

Rc.p - \R~gc.p = 8^ (^V„0V3</) - iv^0V^#„;3 - V{4>)~go.0 + C{4>)f^^ (41) 

This is the first equation for the coupled Einstein-scalar field-matter system. 
Now consider equation ((?T|) and take its trace. We get 



f'{R)R - 2/(i?) + -iUF = 8ttT 



M 
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Multiply through by ■j^J^P-'^ to get 



1 /167r 
16^V ~3~ 



f{R)R-2f{R) 



F2 



If we make the conformal transformation g ^ g and T^^ i-^ as before, we 
find that the wave operator on F transforms like 

UF = FUF - y^FVaF 

Recalling that F = e^^^ then gives 



Then substituting into gH), writing T^^ as FT'^^ (since T^^ f^' = T^ji9"'^ = 
-^^^^Q/i^"'^ ^ F^^T*^), and recalling the definition of the derivative of the 
potential, V{4>), of the scalar field given in ([22]) gives 



which gives equation p5|) . where the coupling constant c((/)) is given by 



1 16n 1 



(43) 



To recover the third (and final) equation, given by p6|l . we must use the 
fact that the energy-momentum tensor of the matter in the f{R) theory is 
divergence-free i.e. 



M 

a/3 







(44) 



The definition of the covariant derivative of a tensor gives 

VrpAI O rj^AI -p(5 rj^AI T^S rpKI 

Again, we can make the conformal transformation g ^-^ g,T^^ i-^ T^^ . We get 
the following 



VarpAI 



5 rfM 



1 



- -F-^d^FS' + d^Fdi - d'F~g^^) T, 



'-50 
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which simphfies to 

and since we require that the divergence of T^p is zero, this gives 

FV^f^l - V^FT^^'p + ]^^pFf^' = 
Since wc know ah-cady that C((/>) ~ F^^, wc sec that 

and so we find that 

Recalhng that F = e^/"^*^ then gives 

which, recalhng the definition of c(0), gives the third and final condition for g 
to be a solution to the Einstcin-scalar field-matter system. 

Therefore, if the metric gap is a solution to the /(i?)-matter field equations, 
we find that under suitable choices of (j) and Vi^cj)) and the coupling constants 
C{(j)) and c(0), given in equations ([M)) . (|1D|) and the metric = 

Fgap is a solution to the coupled Einstein-scalar field-matter system described 
by equations ([M]) . ([55]) and ((55)1 . We can thus apply the results given in [17] 
on the dynamics of solutions to the coupled Einstcin-scalar field-matter system 
and use them to prove results on the dynamics of the /(i?)-matter system. 

4.2 The potential admits a strictly positive minimum 

As in the vacuum case, we will consider first the situation where the potential 
admits a strictly positive minimum. We can use the results of [17] to prove 
results for the /(i?)-matter system. 

We would first like to show, as in the vacuum case, that we can change 
the assumptions used in [17] on the potential of the scalar field to prove an 
analogous theorem for the Einstein-scalar field-matter system. We show that 
the following Theorem holds (analogous to Proposition 3 in [17]): 

Theorem 6. Consider a solution, g, of the coupled Einstein- scalar field-matter 
system, given by equations \S4^ , i35\) . \3b]) . of Bianchi type /- VIII. Suppose that 
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the energy-momentum tensor, Tap, satisfies the assumptions (as in [17]) 

(DEC) jajbg""' < ~P\ Sabr" < 3p 

(SEC) p + Sabg'''' > 

the coupling constants satisfy (as in [17]) 

(C) |c(0)| < CoC{4>) 

and the potential, V , of the scalar field is of class and satisfies ( alternative 
assumptions to [17]) 

(VI ) ¥{(/)) > Vo for a constant Vq > 

(V2) V tends to oo as (p tends to its endpoints. 
If the solution is initially expanding (H > 0) then, as t ^ oo, the following 
limits are attained: 

(i) H ^ Hoc, with Hoo>Hq>Q 

(u) V ^Vo^, where F^o = 

(iii) y ^ 

(iv) 4>,([> Q 

where H is the Hubble parameter for the coupled Einstein- scalar field-matter 
system. 

Proof. First note that the proof of Proposition 3 in [17] first requires proving 
Lemma 1 and Proposition 2 in [17]. To prove these for our alternative assump- 
tions on the potential of the scalar field, we must only find where assumptions 
(PI), (P2) and (P3) were used and check that we get the same results with 
assumptions (VI) and (V2). For Lemma 1, it was only necessary to use (PI) to 
show that V is positive. This is clearly obtained from our alternative assump- 
tion (VI). (P2) was used to give the bound on V'{(f>). From our assumption 
(V2), however, as in the case for vacuum f{R), we see that, since V is bounded 
to the future, we must have that </> remains bounded away from the endpoints of 
the domain oiVa.st~^ oo. Then, since V is smooth, we find that V' must also 
be bounded to the future. The rest of the proof of Lemma 1 goes through as 
in [17]. Proposition 2 docs not require any of the assumptions (PI), (P2); (P3) 
and so the proof of this is identical to that in [17]. With the results of Lemma 
1 and Proposition 2 at our disposal, it only remains to show how Proposition 3 
can be proved with our alternative assumptions on the potential. 

We first note that parts (i) and (ii) are proved in the same way as in [17], 
since they do not require (PI), (P2) or (P3). For part (iii), [17] requires (P3) 
to show that, in the case where converges as t oo, if converges to plus 
or minus infinity, then we still have that V' converges as t — > oo. With our 
assumptions (VI) and (V2), however, we know already that remains bounded 
as t ~* oo, which rules out the possibility that converges to ±oo, and so we 
no longer require this assumption. The rest of the proof then follows as in [17]. 
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It should be noted that the discrepancy in the constants in part (ii) arise from 
the fact that SttG = 1 in [17]. □ 

We are now ready to transfer these results onto the /(i?)-matter system. 
Recall the discussion from the vacuum case in section 3.2, which gave the condi- 
tions on the f{R) theory that are equivalent to the condition that the potential 
of the scalar field blows up at its endpoints. These conditions are labelled (Fl), 
(F2) and (F3). With this in mind, we can now prove the following theorem on 
the dynamics of the /(i?)-matter system. 

Theorem 7. Consider a solution, gap, of the f{R)-matter system of Bianchi 
type I-VIII. Suppose that the energy-momentum tensor, T^^, given by jST] ), 
satisfies 

(1) jaibg"" < P\ Sabg"" < 3p 

(2) p + Sabg"' > 

and that the function f{R) satisfies 

(3) f{R) - Rf{R) < for all f{R) > 0. 

(4) f{R) ~ Rf'{R) satisfies the conditions (Fl), (F2) and (F3) given in 
section 3.2. 

If the solution satisfies H + ^ > initially and exists globally to the future then, 
as t ^ oo, the following limits are attained: 
(i) H ^ Hoc with Hoo> Ho>0 

(n) -^^jf^^^ ^ K where K = 6il^ 
( "V (TUW 

where H is the Hubble parameter for the f{R)-matter system. 

Proof. From the relationships between f{R) and the potential, V{(l)), of 0, us- 
ing the discussion preceding Theorem 2 in Section 3.2, and from the relation- 
ship between the Hubble parameters, we see that, if g is a solution to 
the /(i?)-matter system satisfying assumptions (3) and (4) of Theorem 7, then 
g = Fg is a solution to the coupled Einstein-scalar field-matter system, where 
the potential, V, satisfies (VI), (V2) of Theorem 6. In addition, if the matter 
of the /(i?)-matter system satisfies (1) and (2), then we see from the confor- 
mal transformation described in the previous section that the matter, T, of the 
Einstein-scalar field-matter system satisfies (DEC), (SEC). Lastly, from the def- 
initions of the coupling constants C{(j)) and c{(j)) (given in (03])), we see 
that the ratio between them is constant. Hence the condition (C) of Theorem 
6 is automatically satisfied. 

Therefore, if g is a solution to the /(i?)-mattcr system satisfying (1) - (4) 
of Theorem 7, then g = Fg is a solution to the coupled Einstein-scalar field- 
matter system satisfying (DEC), (SEC), (VI), (V2) and (C) in Theorem 6. So 
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the results of Theorem 6 hold for g. Then parts (ii) and (iii) of Theorem 7 are 
proved directly from the definitions of V{4)) and in terms of i?, f{R) and 

f'{R), although we must still prove the relationship between K and i?oo- 

For part (iii), recall from that H ^ H + ^ and that 4> = \J^^j^- From 
Theorem 6, we know that ^ as t ^ oo. Hence we find that H behaves like 
H at late times and so H ^ Hoo{~ Hoo) with Hoo > Ha > 0. 

To show the relationship between K and Hoo, note that the quantity 

f{R) - Rf'jR) 

is equal to 16ttV and so tends to IdnVcxD which, from Theorem 6, is equal to 
6-^^. Since we have shown above that = Hoo, this completes the proof. □ 

4.3 Intermediate Inflation 

As in the vacuum case, it has also been proven for the coupled Einstein-scalar 
field-matter system [17] that accelerated expansion occurs for the "intermediate 
inflation" case, where the potential of the scalar field is strictly positive, but 
tends to zero as tends to infinity. Once again, we can use results from [17] 
on the dynamics of the Einstein-scalar field-matter system to prove results on 
the dynamics of the /(i?)-matter system. We assume, as in section 3.3 for the 
vacuum intermediate inflation case, that f'{R) — > oo as t — > oo. We prove the 
following theorem: 

Theorem 8. Let g be a solution of the f{R)-matter system of Bianchi type 
I-VHI. Suppose that the energy-momentum tensor, T^^, satisfies 

(1) jajb9'''<p\Sab9'''<ip 

(2) p + SabQ"'' > 

and that the function f{R) satisfies 

(3) f{R) - Rf'{R) < with ^'•^j7(gff^' ^0 as f'{R) oo 

(4) 2f{R) - Rf'iR) < for all R 

(5) -I < limsup < as f'{R) ^ oo 

// the solution satisfies H + j: > initially and exists globally to the future then 
the solution undergoes accelerated expansion at late times. 

Proof. Once again we use the relationship between the /(i?)-matter system and 
the coupled Einstein-scalar field-matter system and the discussion in section 
3.3 of this paper to show that, if 5 is a solution to the /(i?)-matter system 
satisfying the assumptions of the theorem, then g = Fg is a solution to the 
coupled Einstein-scalar field-matter system satisfying the conditions of section 
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4 of [17]. In particular, for the g solution, the result of Proposition 5 in [17] 
proves that, if a := limsup(— V^'/y) < \J~^^- then the following bound is 
obtained: ^ 

The requirement of the bound on the ratio V' /V is deduced from assumption 
(5) of Theorem 8 above. This then gives 

for sufficiently large t which, as we saw in section 3.3, is equivalent to 

H + H'^ >Q 

which gives a > and therefore accelerated expansion at late times. It should 
be noted once more that there is a slight discrepancy between the factors used 
in this paper and those used in [17], due to the notation SttG = 1 in [17]. 

Similarly to section 3.3, we must use the relationship between the Hubble 
parameters of the two systems and information gained during section 4 of [17], 
to show that the expression H+H^ behaves like H + H^ as t ^ oo, which proves 
accelerated expansion for the /(i?)-matter system. This is done by noting that 
the discussion after Proposition 4 in [17] tells us that ^ as t ^ oo. Then 
this, together with the decay results of Proposition 4 and the equation (13) 
of [17] gives a bound on which tells us that — > od as t oo. The bound 
on V' /V and the assumption that y ^ as t — > oo tells us that y — > as 
t — > oo. Then the equation of motion 

4> + 3H<j> + V'((j)) = c{(j))f 

and the information above gives (/) — > as i — > oo. Then (c.f. section 3.3), we 
see that H and H behave like H and as t ^ oo, which gives the result that 
H + > for sufficiently large t, which concludes the proof that the solution 
undergoes accelerated expansion at late times. □ 



5 Conclusions 

In conclusion, we have shown that it is possible to exploit the connection between 
the /(i?) theory and the Einstein-scalar field system to prove mathematically 
rigorous results on the late-time dynamics of the f{R) theory. In particular, 
we have shown that, for a Bianchi type I- VIII space-time, by placing certain 
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restrictions on the function /(i?) of the f{R) theory (obtained according to 
restrictions placed on the potential, V, of the scalar field in the Einstein-scalar 
field system), we have proved that solutions to the f{R) theory satisfying these 
restrictions undergo accelerated expansion at late times. We have shown that 
this is true in both the vacuum case and the case where we add ordinary matter 
to the f{R) theory. 
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